Reconstructing the Peculiar Velocity of the Local Group with Modified Gravity and 

2MASS 



Aidan C CrookQ Alessandra Silvestri, and Phillip Zukirjj] 

Department of Physics, Massachusetts Institute of Technology, 
Cambridge, MA 02139, USA 



O 
O 

(N 
(N 



u: 
43 : 

Oh- 
i < 

o ■ 

c3 . 



(N 

O 

On 
O 



The peculiar velocity of the Local Group, reconstructed from inhomogeneities in the local density 
field, differs in direction and magnitude from the velocity inferred from the Cosmic Microwave 
Background dipole. We investigate whether generalized theories of gravity, which predict a modified 
growth of perturbations, are able to alleviate this discrepancy. We introduce a general formalism for 
calculating the real-space peculiar-velocity field for modified gravity and theories with interactions 
in the dark sector. For different classes of theories — scalar tensor and higher-dimensional gravity 
— we reconstruct the Local Group peculiar velocity using groups of galaxies identified in the 2MASS 
Redshift Survey. We show that, for realistic parameters, modifications to General Relativity cannot 
account for the angular discrepancy between the reconstructed Local Group velocity and the dipole 
in the Cosmic Microwave Background. 



I. INTRODUCTION 

A little over a decade after the discovery of the Cos- 
mic Microwave Background (CMB) [I[, reports emerged 
of a dipole anisotropy in its spatial temperature distri- 
bution 0, 0|. The origin of the dipole anisotropy is 
now widely accepted as a Doppler effect arising from 
the motion of the Sun with respect to the CMB rest 
frame. Measurements by the COsmic Background Ex- 
plorer (COBE) 0, 0], confirmed more recently with 
WMAP, imply a Solar System velocity with respect to 
the CMB of 369.0 ± 0.9 km s" 1 @. Efforts to infer the 
motion of the Sun with respect to the Local Group (here- 
after, LG) centroid [e.g. 0, H[ further imply that the LG 
has a net peculiar velocity in excess of 600 km s^ 1 with 
respect to the rest frame of the CMB, in the direction 
(/,&)- (276°, 30°) [e.g. I EH. 

The large peculiar velocity of the LG is believed to 
have been gravitationally induced. A number of studies 
have attempted to make a comparison between the LG 
peculiar velocity inferred from the dipole in the CMB 
and the velocity reconstructed from inhomogeneities in 
the local density field [e.g., ©, El El GJ, El- The results 
vary significantly depending on both the source of the 
survey data as well as the reconstruction technique used. 

In reconstructing the velocity from local inhomo- 
geneities, one needs to assume a bias factor to infer the 
total mass density from the observed one. For the usual 
assumption of a linear bias, the predicted magnitude of 
the reconstructed velocity depends on the value of the 
bias, while its direction is insensitive to any change in 
the value of the bias; therefore a discrepancy in the angle 
of the reconstructed LG peculiar velocity is of more sig- 
nificance than a discrepancy in the magnitude. The most 
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recent effort by Lavaux et al. [H| predicts a LG velocity 
~ 50° from the CMB direction, with an estimated error of 
22° (95% confidence) from the reconstruction technique. 

There arc a number of possible reasons for this large 
discrepancy. Since the velocity is reconstructed from a 
survey of galaxies, it is natural to wonder whether the 
discrepancy arises because of the finite size of the sur- 
vey. Including the error caused by a finite survey volume, 
as well as the error associated with the reconstruction 
method, the work of Lavaux et al. [l3[ suggests the dipole 
directions should agree to within 25° at 95% confidence. 
Indeed, there are a number of other possible explanations 

from reduced survey completeness in the zone of avoid- 
ance, to invalid assumptions in the estimation of mass, to 
the assumption of linear bias. Unfortunately, due to the 
limitations of the surveys and available data, it is difficult 
to avoid making assumptions of this kind. An analysis of 
each of these postulates is beyond the scope of this pa- 
per. We focus on the assumption that General Relativity 
(GR) is the correct theory of gravity. This assumption 
of GR is common to all previous reconstruction efforts. 

While GR has been tested to high accuracy on small 
scales, such as the solar system [14|, the constraints com- 
ing from cosmological tests are weaker. In this paper we 
explore the effects of large-scale modifications of grav- 
ity on the reconstruction of the LG peculiar velocity 
obtained from an analysis of groups of galaxies in the 
2MASS Redshift Survey [IS 51 E3 El- We consider 
models of modified gravity and interacting dark energy 
that serve as alternatives to ACDM and that modify the 
growth of perturbations in a time- and scale-dependent 
manner. Models of modified gravity that predict a scale- 
independent growth of structure cannot modify the di- 
rection of the LG peculiar velocity with respect to the 
GR prediction, since the peculiar velocity toward each 
attractor would simply be rcscalcd by a common factor 
in this case. However, a scale-dependent modification 
will change the weight of objects at different distances 
from the observer, and is therefore expected to rotate 
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the direction. 

After generalizing the relationship between the pecu- 
liar velocity and density field to allow for modifications 
of GR and interactions in the dark sector, we study two 
classes of models that introduce scale-dependent modifi- 
cations. The first class of models we consider is inspired 
by scalar-tensor, e.g. f(R), as well as coupled dark en- 
ergy theories; generically, they introduce a Yukawa-type 
correction to the Newtonian potential. The second class 
is motivated by extra-dimensional theories 0, H(| and 
introduces a mass scale for the graviton. 

The paper is organized as follows. In Sec. Owe gener- 
alize the relationship between the peculiar-velocity field 
and the density field in real space to models of modified 
gravity and models with interactions in the dark sector. 
We then specialize to the specific cases mentioned above. 
In Sec. IIIII we discuss our reasons for utilizing the 2MASS 
Groups sample and describe the method used to recon- 
struct the LG peculiar velocity from data in the 2MASS 
Redshift Survey. We discuss our findings in Sec. |lVl and 
summarize our conclusions in Sec. IVl 



II. PECULIAR VELOCITIES IN MODIFIED 
GRAVITY 

In the linear regime, assuming zero pressure, the con- 
tinuity equation reads 



V- v 



dS 
i— 
dt 



(1) 



where v is the peculiar velocity field, a is the scale fac- 
tor, 5 = Sp/p is the matter density contrast, p is the 
background density, t is the cosmic time and the spa- 
tial derivatives are taken w.r.t. comoving coordinates. 
For modified gravity, eq. {1} holds generically since the 
energy-momentum conservation equations are left un- 
changed. For models in which there is an additional 
interaction in the dark sector, like coupled quintessence 
models [HI, HH, the conservation equations are modified, 
since the energy-momentum tensor for dark matter is not 
separately conserved. However, in the sub-horizon limit, 
cq. ((T|) holds approximately for dark matter [22j , and we 
can use it for our reconstruction algorithm. 

Assuming vorticity is negligible on the scales we con- 
sider [23j . we can rewrite the peculiar velocity in terms 
of the velocity potential v = V4> v - Plugging this into 
the continuity equation |T]), and transforming to Fourier 
space, we obtain: 



i>v = -grf(k, a)Sk(a) 



(2) 



where H = d\na/dt, k is the comoving wave vector, and 
/(fc,a) = dliid/dhia is the growth factor. 

We can solve for the growth factor by combining the 
Euler equation and the continuity equation ([T]) to find 



the equation for the growth of structure: 
dt V dt 



(3) 



where $ is the (Fourier transformed) time-time compo- 
nent of the perturbed metric in Newtonian gauge. In 
models with a coupling in the dark sector, the r.h.s. of 
cq. ([3]) contains a contribution from the coupling term 
in the Euler equation. As a result, dark matter particles 
respond to an effective potential; however, this effect can 
be absorbed in a modification of the Poisson equation 
which we introduce below. 

We shall now substitute for <t in terms of the density 
contrast. This can be done via the Einstein's equations, 
which in general will be modified. We introduce a generic 
function of time and scale, p k (a), to describe how mod- 
ified theories change the relationship between $ and 5k 
on sub-horizon scales: 



fc 2 $ = 4irGa 2 p(a)p k (a)6 k 



(4) 



Combining equations ^ and and rewriting the 
resulting equation in terms of the growth factor / and 
In a, we obtain a nonlinear differential equation for the 
growth factor: 

df „9 dlnH\ „ 3^ . . , , ,„. 

dh +f H 2+ W / = 2 n - (a)/ife(a) ' (5) 

where fi m (a) = 8irGp m (a)/3H 2 (a). After solving eq. (J5]) 
for /(a, fc), we obtain the velocity potential in real space 
via the convolution theorem: 

d> v (x, a) = a _ [ d 3 x'Sp(a,x')J 7 (x — x! ,a) , (6) 
A-Kp(a) J 

where x is the comoving distance and T is the inverse 
Fourier transform of 4-7r/(fc, a)/k 2 . Finally, by taking the 
gradient of ([6|) w.r.t. x' we obtain the velocity in real 
space: 



v(x, a) 



aH 
4irp(a) 



d x'Sp(a, x')J- v (x — x' , a), (7) 



where T v = V.F. 

In our reconstruction algorithm we first parametrize 
different theories via the function pk{d), then solve 
cq. © for f(k,a = 1) subject to the initial condition 
/(fc, 10~ 3 ) = 1, and perform the inverse Fourier trans- 
form of f(k,a = l)/k 2 to obtain T(x — x',a = 1) and 
T v . The initial condition for the growth factor is taken 
to correspond to the GR solution during matter domi- 
nation. In other words, the theory must reduce to GR 
at early times. The reconstruction methods based on 
GR commonly use the approximation f(a) = r2 m (a) 6 / 11 
[e.g. |24|; in our analysis we compare the predictions of 
modified models of gravity to the GR result obtained by 
setting /ifc(a) = 1. 
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While the equation for the growth factor ([5]) can be 
solved numerically, we can also proceed analytically with 
some approximations. Although the analytical method 
has limited applicability, as we will discuss, it is very 
useful in understanding the effects of modified gravity. 
Moreover it works well in the relevant range of param- 
eters. In the remainder of this section we present the 
analytical results and then describe the two classes of 
modified gravity models that we consider in our analysis. 

The equation for the growth factor can be solved by 
quadratures if one neglects terms of 0((f — l) 2 ). The 
approximate solution is given by [25|: 

f(k, a)~J^J Q da ' a' z H(a') + ^Q m (a')fi k (a') 

'(8) 

Expression ([8]) is a good approximation when / does 
not depart significantly from unity throughout the expan- 
sion history. It systematically overestimates the actual 
growth factor. For ACDM (fi k = 1), taking Q° n = 0.3, 
it overestimates the growth factor today by 7%. For the 
models of modified gravity that we consider, the preci- 
sion of © depends on the scale and on the parameters 
of the theory. This will be discussed in the following 
subsections, after each model is introduced. 

Given the above analytical approximation for the 
growth factor, we inverse Fourier transform Airf/k 2 to 
obtain the following approximate Green's function for the 
velocity potential. 



T{x-x',a) ps — L f da'a' 3 H(a') 
a H J 



1 



\x — x 



+-n m (a')g(x-x' > a') 



(9) 



In the above, Q{x — x', a') is the inverse Fourier trans- 
form of 4iriik(a)/k 2 . Finally, the approximate Green's 
function for the peculiar velocity is: 



J- V (x — x , a) 



a 4 H 



da'a' 3 H(a') 



x — x 

ry /-*»/ 1 3 



+ -0, m {a')Q v {x - x', a') 



(10) 



where Q v = V<7 . It is also informative to solve cq. ^ for 
the exact Newtonian potential. One finds: 



a) = -Ga 2 / d 3 x'Sp(x')g(x - x', a) . (11) 



Eq. (|T0|) allows us to understand features of the nu- 
merically reconstructed LG velocity dipole for gravity 
models that introduce time- and scale-dependent mod- 
ifications to the growth of structure. We focus on mod- 
els that tune their expansion history to be indistinguish- 
able from ACDM. Therefore, for all the models we con- 
sider, H will be parametrized by the ACDM solution with 



the 5-year WMAP results for the cosmological parame- 
ters, i.e. fl^n = 0.27, H = 70.5 km s -1 Mpc _1 , and 
fl A = 0.73 [2ff . 

In the remainder of this section we specialize to two 
broad classes of models which encompass alternatives to 
the ACDM scenario in the context of cosmic acceleration. 
The first class is scalar-tensor theories and dark energy 
models with interactions in the dark sector. The second 
class is extra dimensional theories of gravity. 



A. Scalar-tensor theories and Yukawa-type 
interactions 

To describe scalar-tensor theories, e.g. f(R) the- 
ories, and models of coupled dark energy, we con- 
sider the parametrization introduced by Bertschinger and 
Zukin [27]. The rescaling of Newton's constant is of the 
form: 



Mfc(a) 



1 + aifcV 
1 + aik 2 a s 



(12) 



In the above, 011,0.2, and s are free parameters. The 
parameter s is dimcnsionlcss while a\ and oti have units 
of length squared. In order to recover the GR solution at 
early times, we take s > 0. We take a.2 > so that the 
coupling is finite, and we take a\ > so that gravity is 
attractive. In the specific case of f(R) theories, the ratio 
ai/«2 = 4/3 and s <~ 4. 

Defining £ = x' — x and £ = the Green functions 
are: 



G&a) 



1 

I 

e 



a i 

(12 



(13) 



(14) 



The above expression for Q, together with eq. (fTTj) . re- 
veals that the Newtonian potential, at late times, con- 
tains an additional Yukawa coupling with (comoving) 
mass m e ff = l/y / a2d s . The Yukawa correction is ex- 
pected since alternative theories typically introduce ex- 
tra degrees of feedom which mediate a fifth force between 
DM particles. In the limit that £ m e g- ^> 1, the Yukawa 
coupling is suppressed and we recover GR. The effective 
coupling, ol\ /a.2 — 1, is positive when the fifth force me- 
diator is a scalar field, which is the case for scalar-tensor 
theories and coupled quintessence. However, it can be 
negative if the force is mediated by a vector field or a 
phantom scalar field [2^|. Astrophysical and cosmologi- 
cal measurements set combined constraints on the length- 
scale and coupling of the fifth force; for scales 100 kpc < 
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< 10 Mpc, the intracluster gas distribution [29| and 
cosmological data [30[ give 0.5<ai/a2<2.3. 

We quantify the error of our analytical results by com- 
paring the approximate Green's function (eq. I10[) with 
the exact Green's function defined in eq. ([7]). The error 
depends on the model parameters and the scale £; it can 
be traced to how closely / ~ 1 throughout the expansion 
history. On large scales, [i^ ~ 1 over most of the expan- 
sion history, and the analytical Green's functions over- 
estimate the exact solution by ~ 7%. On small scales, 
the errors are smaller (larger) than 7% for a±/a2 > 1 
(< 1) since the growth factor is larger (smaller) than the 
GR solution. For ai/a 2 < 1, the error becomes larger 
on small scales for larger y/ai and smaller s, since these 
all contribute to smaller growth factors. The error for 
y/ai = 100 Mpc, ai/o-2 = 0.5, and s = 2 at £ =10 Mpc 
is ~ 30% — a factor of 3 larger than typical errors on 
mass and distance estimates. Hence, the analytical re- 
sults are accurate to ~ 7% on scales 5 < £ < 150 Mpc 
for ot\/ot2 > 1) but should not be trusted on small scales 
for u\/a.2 < !• 



B. Higher-dimensional gravity 

Next we consider higher-dimensional extensions of the 
Dvali-Gabadadze-Porrati fDGP) I31J scenario, within the 
cascading gravity framework of \19L l20j. In this scenario, 
the matter fields are constrained to a 4-dimensional brane 
embedded in a higher dimensional bulk of infinite volume, 
into which the graviton can leak. The 4D gravity carri- 
ers are massive gravitons, with a longitudinal mode that 
mediates a fifth force. The grow th of structure for these 
models has been studied in [32l, [33j ] , where it was found 
that the Poisson equation is modified by the time- and 
scale-dependent factor 



We specialize to the case D = 6, which implies a = 
[H. We find: 



[J-k(a) 



I - g(a) 



1 + (,002(1- 



- kr c • 



-a) 



(15) 



where 



sO) = - 



D-A 



D 



2 J l + 2(ffr c ) 2 ( 1 -«)(l + ^) 



(16) 



D is the number of spacctime dimensions, < a < 1 and 
depends on the number of dimensions, a prime denotes 
derivates w.r.t. In a. and r c is the physical mass scale 
of the graviton. In [32|], it was found that the relevant 
range for the mass scale is r c = 300 — 600 Mpc. Expres- 
sion (TT5|) describes the combined effects of the massive 
graviton and the longitudinal mode of the theory. At late 
times (Hr c < 1) and on sub-horizon scales, the growth 
is enhanced on scales r < r c , by the fifth force mediated 
by the longitudinal mode. On scales r > r c the growth 
is suppressed by the massive graviton. 



S(t,a) = ^(l-g(a))e- a ^, 
Sv&a) = i(l- ff (a))(l + £)e-^ 



(17) 
(18) 



The above expression for (?, together with eq. pip , re- 
veals that the Newtonian potential is modified in ampli- 
tude by (1 — g (a)) and is exponentially cut off on large 
scales. The effective (comoving) mass for the coupling 
is given by m e ff = a/r c . For £ m e ff 3> 1, the Newtonian 
potential exponentially decays. This behavior is different 
from that scalar tensor theories, where in this limit the 
l/£ behavior of GR is recovered on larger scales. This ex- 
ponential decay decreases the range over which gravity 
acts. 

As mentioned above, the error in the Green's function, 
depends on the model parameters and the scale £. The 
analytical approximation works very well in the relevant 
range r c = 300 — 600 Mpc. The error in the Green's 
function for r c = 300 Mpc is 3% at £ = 5 Mpc and 7% 
at £ = 140 Mpc. The error gets larger in magnitude for 
smaller r c ; the analytical results should not be trusted 
for r c < 50 Mpc (where the error at £ = 150 Mpc is 
- 15%). 



C. Summary 

As expected, each of the models introduces scale de- 
pendent corrections to the Newtonian potential and the 
peculiar velocity field. From the Green's functions de- 
rived analytically, it is clear how a scale-dependent rescal- 
ing of Newton's constant, /Xfc(a), modifies the spatial de- 
pendence of the peculiar velocity Green's function. For a 
given choice of parameters, the contribution of masses 
at different distances will be modified in a distance- 
dependent way. As a consequence, the direction of the 
reconstructed LG dipole will be rotated with respect to 
the GR result. The amount and direction of this change 
will depend on the parameters of the models as well as the 
distribution of masses in the survey. Therefore, it is not 
straightforward to estimate the amplitude of the change 
in direction from simple arguments applied to eq. (|10[) . 



III. RECONSTRUCTING THE LOCAL GROUP 
PECULIAR VELOCITY FROM 2MRS 

In this section we first provide an overview of the data 
to be used in the reconstruction (Sec. IIII Al and IIII Bp . 
then discuss the method used to compute the peculiar 
velocity of the LG (Sec. InTC) . 



5 



A. 2MRS Groups Catalog 



The Two Micron All-Sky Survey [2MASS, uni- 
formly mapped the entire sky in the J (1.25 /xm), H 
(1.65 /im), and K s (2.16 pm) bands. By observing in 
the near-infrared, 2MASS is less susceptible to the se- 
vere extinction at low galactic latitudes from which ob- 
servations at optical wavelengths suffer. Since galaxies 
spectra peak at ~ 1.6^m, 2MASS maximizes the num- 
ber of galaxies detected at a specified flux limit, resulting 
in the most complete all-sky survey per formed to date. 
ZRj45-selected samples [e.g., dH, [36l |37|| also avoid the 
extinction problems at low galactic latitudes; however, 
since these samples are based on fluxes in the far-infrared, 
they are biased toward star-forming galaxies and under- 
sample early-type galaxies; we therefore expect models 
derived from such samples to underestimate the masses 
of large clusters. 2MASS makes every effort to obtain an 
unbiased, uniform sample of extended sources across the 
entire sky; combined with the ability to detect galaxies 
at low galactic latitudes, this makes 2MASS an ideal re- 
source to understand the flows in the nearby Universe. 
In this paper we focus purely on the motion of the LG. 

The follow-up 2MASS Redshift Survey [2MRS,E0J] 
is 99.9% complete to K s < 11.25, b > 10°, and provides 
redshifts for galaxies selected from the 2MASS extended 
source catalog. Efforts to reconstruct the LG dipole from 
galaxies in 2MASS [| and 2MRS [HO have produced 
results that are both inconsistent with the CMB and de- 
pend on the reconstruction method. 

The relatively recent work by Erdogdu et al. [l(| at- 
tempts four different reconstructions for the LG dipole. 
While the flux-weighted reconstruction in the LG frame 
suggests a LG dipole 21° ±8° from the CMB result, when 
performing the reconstruction in the CMB frame, the 
result increases to 26° ± 8°. Flux- weighted reconstruc- 
tions necessarily assume an intrinsic mass-to-light ratio, 
and therefore do not account for variation in dark-matter 
bias in clusters compared with isolated galaxies. The re- 
constructions using a number-weighted technique in the 
same paper suggest a LG velocity more than 35° from the 
CMB, and since the majority of distances are estimated 
from redshifts, the appropriate choice for the reference 
frame becomes unclear. The more recent reconstruction 
from 2MRS by Lavaux et al. [H[ predicts the LG velocity 
~ 50° from the CMB direction at lOO/i" 1 Mpc, with an 
estimated error of 22° (95% confidence) from the recon- 
struction technique. It is evident that the assumptions 
built into the method chosen to reconstruct the LG ve- 
locity have a noticeable impact on its direction. 

In addition to errors associated with reconstruction 
methods, there are also errors caused by the finite-size of 
the survey. Assuming a ACDM cosmology, the amplitude 
of peculiar velocity fluctuations for a sample complete to 



140 Mpc 1 is approximately 60 km s _1 . For a ler fluctua- 
tion, given a velocity with respect to the CMB rest frame 
of ~ 600 km s _1 , this implies at most a 6° discrepancy 
for a survey limited to the same radius. Therefore the 
power on large scales neglected due to the finite survey 
volume is unlikely to be the reason for the discrepancy 
in direction. 

Since the distribution of dark matter drives the inho- 
mogeneities in the density field, we want to use the best 
available tracers of dark matter to infer mass. While as- 
suming that baryonic matter traces dark matter provides 
a reasonable starting point, to further assume that the 
mass of baryonic matter is proportional to the mass of 
dark matter ignores any differences between dark mat- 
ter in isolated galaxies or small groups and massive 
galaxy clusters. To avoid making this assumption, we 
will use dynamical mass estimates obtained from galaxy 
groups: we use the catalog of galaxy groups derived from 
2MRS [13, EH as the basis for our reconstruction. Two 
catalogs of groups, constructed using a variablc-linking- 
lcngth percolation algorithm, are presented in the afore- 
mentioned work; in this paper, we use groups from the 
revised 2MRS HDC catalog (constructed with a mini- 
mum density contrast Sp/p = 80); the reason for this 
choice is that these groups are less prone to interloper 
contamination than the alternative (LDC) catalog (e.g., 
[HI), and we therefore expect the dynamical mass esti- 
mates to be more robust. This choice, however, comes 
at the expense of missing some member-galaxies. The 
catalog of groups is truncated at a distance of 140 Mpc, 
where the number of galaxies in the 2MRS K s < 11.25 
catalog between radii r and r + dr falls to half its max- 
imum value. The method we use to reconstruct the LG 
peculiar velocity is discussed below. 

In this work we make an effort to reduce and optimize 
the assumptions in the reconstruction method. However, 
the purpose of this work is to understand the effect of 
modified gravity on the reconstructed LG peculiar veloc- 
ity. And, while choosing a different method, with differ- 
ent assumptions, will undoubtedly change the direction 
of the reconstructed dipole, the conclusions presented in 
this paper are robust, and can be reproduced using the al- 
ternative number-weighted reconstruction technique uti- 
lized by [nj. 



B. Modifications to the catalog 

Since no algorithm to identify groups using only 3 com- 
ponents of phase-space information can guarantee finding 
virialized groups, we have to understand the limitations 
of the data. While the algorithm used to construct the 
2MRS HDC catalog has been tested and optimized using 



1 We choose 140 Mpc since this is the size of the sample we will 
study in this work. 
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simulations (e.g, [39|, 50[)> the groups are still susceptible 
to interloper contamination. 

By examining the sky-distributions and velocity dis- 
persions of the members of the groups with the largest 
contributions to the LG acceleration, we make the follow- 
ing modifications to the group catalog before continuing: 

1. A group behind the Virgo cluster at 
(a=12 h 22 m , 5=5°) with a LG-centric redshift 
of ~ 2200 km s _1 has been merged with the 
Virgo cluster; we manually separate the group 
by requiring that galaxies assigned to the Virgo 
cluster inside an ellipse (a = 3.5°, 6 = 2°, with the 
semimajor axis inclined 30° from North toward 
West) with redshifts greater than 1850 km s™ 1 arc 
re-assigned to a distinct group. 

2. In the group catalog, the nearby M94 group (z ~ 
300 km s _1 ) has been merged with both the M106 
group (z ~ 500 km s _1 ) and the Ursa Major group 
(z ~ 1000 km s _1 ). We manually separate this 
group as follows: We consider those galaxies in- 
side an ellipse (a = 15°, b = 10°, with the semi- 
major axis inclined 60° from North toward East) 
centered on (a=12 /l 40 m , (5=41°); those with he- 
liocentric redshifts below 400 km s _1 are assigned 
to the M94 group, and those with redshifts be- 
tween 400 and 700 km s _1 are assigned to the Ml 06 
group. The remaining galaxies now appear visu- 
ally as 3 distinct groups: those inside the ellipse 
(a = 9°, b = 6°, with the semimajor axis in- 
clined 70° from North toward West) centered on 
(a=12 /l 24 m , (5=29°) are split into two groups, with 
heliocentric redshifts above and below 900 km s -1 . 
Those outside the ellipse comprise the Ursa Major 
group. 

The estimated distances to and masses of the groups are 
then recomputed using the updated list of members. 



Method 



smaller groups due to the increased errors in their mass 
estimates. 

We allow for a linear bias, b, between the total matter 
density and the density in galaxy groups, i.e. 6p g / p g = 
bdp/p, where the subscript g denotes galaxy groups. We 
further assume that mass concentrations are point-like 
and located at the distances specified in the catalog: 
Sp g = MSoir) — p g - Therefore from eq.® the resulting 
velocity vector of the LG (which we define to be located 
at the origin) today is given by 



vlg = 



E 



iirpgb 



1) 



(19) 



where M? ff is the effective mass (see below) of the i th 
group, located at distance \ fl\ 3 and T v is the numerically 
computed Green's function defined below eq. [7] The nor- 
malization of eq. (|19|) depends on the mean matter den- 
sity in galaxy groups, p g , as well as the bias parameter, 
b. For the remainder of this work, we will set 6=1, since 
this is just a scaling parameter that will modify the mag- 
nitude of the LG peculiar velocity by a constant factor, 
but will not affect the direction. 

Since we are using data drawn from a magnitude- 
limited survey, the data are subject to the selection bias 
of the sample (i.e. fewer galaxies are observed per unit 
volume at larger redshifts). Since the sample is restricted 
to small redshifts, we expect that the intrinsic mean mass 
density does not vary with position, however the mass 
density we infer by estimating the masses of observed 
groups will inevitably decrease with redshift. In an effort 
to reduce the impact of the selection bias on our results, 
we apply a correction to the estimated group masses such 
that the effective mass per unit volume does not change 
with redshift. 

We assume the observed mass density has the func- 
tional form 



Pg( r ) = Pg 



■kr 



(r/ro) 2 



J x a+1 e x dx 



(20) 



The 2MRS HDC group catalog (see Sec. HEl above) 
is a powerful tool for reconstructing the LG peculiar ve- 
locity. However, since we are dealing with a catalog con- 
structed from observational data, we must be careful to 
avoid its caveats. The catalog was constructed using a 
variable linking-lcngth percolation algorithm, based on 
that of [4l[, in an effort to prevent an artificial reduc- 
tion in the number of groups identified at large distances. 
The resulting selection function 2 of the group catalog will 
therefore differ from that of 2MRS. In our analysis, we 
consider only groups with 5 or more members, neglecting 



2 Wc define the selection function as the number of galaxies /groups 
observed per unit volume as a function of redshift. 



where a and tq are free parameters. A function of this 
form is motivated by assuming the luminosity density is 
independent of spatial position, and the luminosity func- 
tion has the form of a Schcchtcr function [IH; if the 
mass density is proportional to the luminosity density, 
then the observed mass density would follow the form of 
eq. (f20|) . While wc do not require that these assumptions 
hold true, the functional form satisfies the properties that 



3 Since the groups under consideration are at redshifts below ~ 
0.03, the difference between physical and comoving distances is 
smaller than the errors on the estimated distances to groups. As 
such, we do not apply corrections to the distances and masses 
reported in |l7l Il8l | to account for the cosmology and variation 
in scale-factor at the different groups. 
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FIG. 1: Mass in observed groups. The solid line shows the 
mass (as a fraction of the total mass) , estimated dynamically, 
in observed groups of galaxies with 5 or more members inside 
a sphere of radius r, as a function of r. The dotted line shows 
the function of form (|2ip obtained with a least-squares fit. 
The deviations from the fit for nearby distances are due to 
the overdensity in the Local Supercluster. 



the observed mass density (i) remains finite at the ori- 
gin, and (ii) decays at large distances, and therefore it 
provides a good empirical fit to the data. 

The observed mass inside a sphere of radius r is given 

by 

M{r) = 4tt / p g (r')dr' (21) 
Jo 

and p g is specified by requiring 

M(D) = M tot 

where M to t is the total mass in observed groups, and 
D is the extent of the sampled volume. We fit eq. (|2Tj) 
to M(r) constructed from the catalog of galaxy groups 
(see Figured]), and obtain a = 0.47, ro = 65 Mpc. The 
corresponding p g = 1.24 x 10 10 M Q Mpc -3 . By assigning 
an effective mass to each group, given by 



M 



eflf 



Pg{ r i)/Pg 



(22) 



where the Mf are the projected masses [43[ of the groups 
taken from the catalog, 4 we are assuming that mass we 
do not observe is correlated with the observed galaxy 
groups. 



4 The projected mass estimator assumes that GR is valid on scales 
of order < 0(l)Mpc. The models of modified gravity we con- 
sider, display a non-linear mechanism which restores GR on small 
scales (for which Sp/p > 1), therefore it is reasonable to assume 
standard GR in the dynamical estimates of masses. 



IV. RESULTS 

Before we discuss the effects of modified gravity on 
the direction of the dipole, it is helpful to analyze the 
structure of the reconstructed dipole under GR. The con- 
tributions from nearby structures can be understood by 
reconstructing the dipole using only the mass inside a 
sphere of radius r. The magnitude of the peculiar veloc- 
ity plotted as a function of the sphere radius, is shown by 
the solid line in Figures [4] and [71 As expected, the form of 
the galactic cartesian components expresses the same fea- 
tures as the reconstruction from the 2MRS galaxy sam- 
ple [m inside 10,000 km s 1 . However, since the galaxies 
have been grouped, we observe a discretized version of the 
galaxy-based reconstruction, with slight modifications to 
the amplitude because we use dynamical mass estimates. 

The key features to note are the following. The M94 
group is the nearest significant group, at a distance of 
~ 4 Mpc located in Ursa Major at (l,b) w (149°, 77°); 
lying near the galactic north pole, this group contributes 
significantly to the z-component of the LG dipole. At 
~ 15 Mpc, the Ursa Major group (l,b) « (145°, 67°) pro- 
vides a further attraction toward the galactic pole. The 
Virgo cluster (2MRS HDC Group #720) at (/,&) w (285°, 
73°) at a distance of ~ 16 Mpc provides the largest single 
contribution to the LG velocity of any group, increasing 
the LG peculiar velocity in the z-direction and decreas- 
ing the y-component. Centaurus (2MRS HDC Group 
#730) at a distance of 57 Mpc at (/,&) « (302°, 22°) re- 
duces the y-component and increases the x— and z— com- 
ponents. The Perseus-Pisces cluster (2MRS HDC #219) 
at a distance of 77 Mpc at (l,b) m (150°, -13°) and the 
Norma Cluster (A3627, 2MRS HDC #941) at 72 Mpc, 
(l,b) « (325°, -7°) both drive down the z-component of 
the LG peculiar velocity, but the "tug of war" in the x- 
and ^-directions is clearly evident as the pull of Norma 
is quickly countered by the opposing Perseus-Pisces. 

At larger distances, the volume increases more rapidly 
with distance and the discrete changes in velocity be- 
come less evident. The Coma Cluster (A1656, 2MRS 
HDC #745) at 106 Mpc, (l,b) w (59°, 88°) provides the 
largest single contribution beyond Perseus-Pisces, 5 but 
the dipole appears to have converged by 120 Mpc, some 
46° from the CMB dipole. To assess the propagation 
of errors from the mass and distance estimates we use a 
Monte Carlo technique, assuming Gaussian errors of 10% 
on the distance to the group, in addition to a Poisson er- 
ror for the dynamical mass estimate; we find the error 
on the reconstructed angle is 7° (10°) at 68% (90%) con- 
fidence. Furthermore, A-body simulations performed to 
assess the applicability of linear perturbation theory, sug- 
gest that the uncertainty in the reconstruction method 
is w 16° (23°) at 68% (90%) confidence 0. Combin- 



Note that the Shapley Supercluster at 
sampled volume. 



180 Mpc is outside the 
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ing these error estimates, the discrepancy with the CMB 
dipole should not exceed 25° at 90% confidence, and the 
probability of obtaining an offset of 46° or higher is less 
than 3%. 

It is important to be aware that the conclusions pre- 
sented below will not change if constructed using an al- 
ternative source catalog (e.g. galaxy-based sample) since 
the general form of the contributions toward the LG pe- 
culiar velocity remains the same. 

Armed with this insight, we now turn to discuss the 
effect of modified gravity on the LG peculiar velocity. 
With a LG dipole under GR directed > 40° away from 
the CMB expectation, it is clear that modified gravity 
must be able to rotate the direction through an angle 
of this order of magnitude to provide a solution to the 
problem. 

A. Scalar-tensor theories and Yukawa-type 
interactions 

Using the scalar-tensor paramctrization shown in 
eq. (Tj"2"]) . we numerically calculate the LG peculiar ve- 
locity vector for theories that introduce an additional 
Yukawa interaction. We compute this vector for s = [2, 4] 
and ai/«2 = [0.5, 1.5, 2.3], as a function of ^/a^. The 
values of ot\l u-i are chosen in order to be consistent with 
current limits [2£|[3(|, as discussed in Sec. [Til This range 
of s is expected for typical Chameleon-type models |45| , 
and the range of yfo2 is chosen so that GR is recovered 
at early times (i.e. fXk — ► 1). In FigurcEJ we plot the an- 
gular deviation of the reconstructed LG velocity from the 
GR result for this class of models; the amplitude of the 
reconstructed velocity, relative to the GR case, is plot- 
ted in Figure [3] From Figure [21 it is evident that the 
maximal angular deviation is ^3°, which is significantly 
smaller than the ~ 40° needed to align the LG velocity 
with the direction inferred from the CMB dipole. 

The parameters s, ai/a2, and «2 have different ef- 
fects on the velocity dipole. The parameter s changes 
the time at which modified gravity becomes dominant. 
For larger values of s, modified gravity becomes relevant 
at later times, leading to smaller deviations from GR. 
Larger s also shifts the features in Figures [5] and [3] to 
larger values of ^Jai since a Yukawa coupling which acts 
over a larger range compensates for the smaller time du- 
ration over which modified gravity is relevant. In order 
to understand the effects of oli and 0.1/0.2, let us look 
at the Green's function (|14[) derived analytically. As dis- 
cussed previously, the approximation we have used in de- 
riving (|14[) works well for values of ai/<X2 > 1 (and on 
large scales for ax/a2 < 1) and is useful for interpreting 
the numerical results. 

Consider the deviation caused by modified gravity de- 
fined as the fractional difference in the magnitude of the 
velocity, 



"GR 



and assume s is fixed. Looking at the approximate 
Green's function (eq. [TJJ, we see that m e g sets the tran- 
sition scale between two asymptotic behaviors: Av oc 
{pt.i/ot.2 — 1) for £m c ff <C 1 and Av ~ for £m c g 
1. Hence, calculating the LG peculiar velocity effec- 
tively becomes a weighted average, with the contribu- 
tion of a mass to the LG peculiar velocity weighted by 
~ 0.44(ai/a2 - 1) + 1 for £,m ctf < 1, and by — 1 for 
£m e ff 3> 1; both weights are relative to the predicted 
GR contribution. Since m e g = 1/ \/a2a s , changing 02 
changes the scale at which the shift in weighting occurs 
and changing a.i/a.2 changes the weighting itself. 

Note that the sign of (ai/tt2 — 1) determines whether 
the contribution to the LG peculiar velocity from galaxy 
groups on scales £?n G g- < 1 is smaller or larger than the 
GR case. In the case of ai/«2 > 1, the change in di- 
rection brings the LG dipole closer to the CMB-implied 
direction, whereas in the case of a±/oi2 < 1, the shift 
is away from the CMB dipole. Since the amplitude of 
the change in angle is significantly less than the error 
on the reconstruction, we cannot draw any conclusions 
from this observation. Since our sample probes 5 Mpc 
< £ < 140 Mpc, the above analytically computed limits 
for Av are roughly reproduced in Figure[3]for ai/012 > 1. 
For ai/aa < 1, the analytic limit for £m e g <C 1 does 
not match the numerical result since the approximation 
breaks down. 

Figure [2] shows that the maximal angular deviation 
from the GR reconstructed dipole occurs for in the 
range 20-40 Mpc (for s = 2) and 30-50 Mpc (for s = 4), 
depending on the choice of a.\Ja.2- This maximum devia- 
tion is ~2-3°, significantly less than the ~40° required to 
explain the discrepancy with the CMB. The angular de- 
viation goes to zero at the endpoints of the domain since 
the dipole direction is insensitive to a constant rescaling 
of all velocities. The value of m g at which Figure [2] 
peaks, however, depends on the distribution of galaxies. 
This is analyzed below by reconstructing the LG pecu- 
liar velocity using only groups inside concentric spheres 
of increasing radii, which is shown in Figure [H This al- 
lows us to identify which groups drive the amplitude and 
direction of the peculiar velocity, and how this changes 
under specific cases of the Yukawa- type. In what follows, 
we specialize to the parameter values used for Figure 0] 
and discuss the details of the plots. 

The left panel shows the LG peculiar velocity for 
011/012 = 0.5 (repulsive Yukawa), s = 2, and y/a2 = 30 
and 100 Mpc. For ^/a2 = 30 Mpc, the contribution 
to the peculiar velocity from the groups inside 30 Mpc 
(dominated by the M94 and Ursa Major groups, and the 
Virgo Cluster) is significantly reduced compared to the 
GR case (since oti/a.2 < 1), and the influence of groups 
beyond 30 Mpc (most notably Centaurus, Perseus-Pisces 
and Norma) are dampened when compared with GR. As 
such, the relative contribution of Centaurus is less than 
in GR, so we don't see the negative shift in y or positive 
shift in x at a distance ~ 55 Mpc. At very large dis- 
tances, the increased weighting of the groups, and larger 
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FIG. 2: Angular deviation of the LG peculiar velocity reconstructed under Yukawa-type modified gravity (or dark- interactions), 
from the direction of the standard GR reconstruction. The left panel shows the case s = 2 and the right panel s = 4, for three 
different values of the effective coupling ai/aa: 0.5 (solid line), 1.5 (dotted line), and 2.3 (dashed line). 




FIG. 3: Deviation in magnitude of the reconstructed LG peculiar velocity relative to the GR result, under Yukawa- type modified 
gravity (or dark-interactions) . The figures show Av (eq. I23[l as a function of 1/02 . The left panel shows the case s — 2 and the 
right panel s = 4, for three different values of the effective coupling 0:1/02: 0.5 (solid line), 1.5 (dotted line), and 2.3 (dashed 
line) . 



volume of the concentric shells, means the contribution 
to the peculiar velocity from masses within the shell is 
very sensitive to the distribution of mass, and a depar- 
ture from GR is expected. However, the amplitude of the 
contribution to the peculiar velocity from these shells is 
not sufficient to significantly alter the direction of the 
dipolc. Since the inner groups dominate the contribu- 
tion toward the galactic pole, the net result is a small 
shift in the direction toward the galactic plane, bringing 
the LG dipole ~ 2° farther from the direction expected 
from the CMB. For y/ct2 = 100 Mpc, the transition scale 
is close to the limit of the survey. The resulting peculiar 
velocity therefore follows the GR result, reduced in am- 
plitude by a constant factor, up to the transition scale. 
Since the dipolc has converged by ~ 120 Mpc, there is lit- 
tle angular difference from the GR result. Hence, while 
the amplitude of the LG velocity is reduced, the angle 
remains consistent with GR to ~ 1°. 



The right panel shows the case of 011/0:2 = 2.3 (attrac- 
tive Yukawa), s = 2, and = 30 and 100 Mpc. We see 
the reverse of the above scenario. Since a\/a2 > 1, the 
contribution to the peculiar velocity of groups inside y/a2 
is amplified. Hence for 1/02 = 30 Mpc, the increased con- 
tributions from the nearby M94 and Ursa Major groups 
and the Virgo Cluster increase the z-component of the 
velocity over the GR result. For v / cv2 = 100 Mpc, the 
contribution from all groups is amplified and the graph 
follows the GR result up to ~ 100 Mpc, after which 
the relative contributions of more distant groups is sup- 
pressed, most significantly - the Coma Cluster - which 
dampens the jump in positive z. The groups beyond 100 
Mpc have large effective masses (due to the correction 
for the selection bias) , and one would expect that chang- 
ing the weighting of these groups via modified gravity, 
would have significant impact on the reconstructed ve- 
locity. However this effect competes with the isotropy 
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FIG. 4: Components of the LG peculiar velocity as a function of the size of sphere-of-infiuence under Yukawa-type modified 
gravity. The graphs are constructed using s = 2. The left panel shows the case a.x/a.2 = 0.5 and the right panel 01/02 = 2.3. 
The graphs show the cartesian components of the LG peculiar velocity in galactic coordinates. The vertical axis is normalized 
such that at the limiting distance of the sample (140 Mpc), the magnitude of the GR-case velocity is unity. The solid line 
represents the GR reconstruction; the dotted and dashed lines represent the peculiar velocities reconstructed for two different 
scale lengths: ^fai = 30 Mpc and 100 Mpc, respectively. Note that the vertical axes of the different cartesian components do 
not have the same scales. 



on these scales; indeed, changing the weighting of an 
isotropic distribution cannot alter the direction of the 
dipole. Hence the effect of modified gravity on the contri- 
bution toward the LG peculiar velocity from the distant 
groups is somewhat suppressed. 



B. Higher-dimensional gravity 

Using the parametrization shown in eq. (|18[) . we nu- 
merically reconstruct the LG peculiar velocity for extra- 
dimensional theories. Since we have chosen D = 6, this 
leaves the crossover scale r c as the only free parameter. 
While the relevant range is 300 Mpc < r c < 600 Mpc p^ . 
we show the results for 10 Mpc < r c < 600 Mpc so the 
effect of this formulation is evident to the reader. Figure 
[5] shows the angular deviation of the reconstructed veloc- 



ity from the GR result, plotted as a function of r c ; the 
fractional difference in the magnitude of the velocity is 
shown in Figure [SJ 

The parameter r c causes deviations from GR in two 
distinct ways. To start with, it sets a transition time 
(Hr c < 1) after which modifications to GR become im- 
portant. As equation (fl6|) shows, g(a) evolves from ap- 
proximately zero at early times, to a value of ~ — (D — 
4)/(D — 2) at late times, (when Hr c <C 1). Larger 
r c causes modified gravity to become relevant at later 
times, leading to smaller deviations from GR. Moreover, 
although we restrict ourselves to D = 6, it is clear that 
more dimensions will cause larger deviations from GR. 

The parameter r c also sets a transition scale. At late 
times (Hr c <C 1), and on sub-horizon scales, the growth 
of structure is enhanced for k > a/r c , while it is sup- 
pressed for k < a/r c . This behavior is reflected in the 
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FIG. 5: Angular deviation of the LG peculiar velocity re- 
constructed under degravitation, from the direction of the 
standard GR reconstruction 
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FIG. 6: Deviation in magnitude of the reconstructed LG pe- 
culiar velocity for the higher-dimensional gravity case, relative 
to the GR result. The figures show Av (eg. I23p as a function 
of r c . 

approximate Green's function (|18[) . where contributions 
from masses at £ S> r c /a are suppressed. Masses at 
£ "C r c /a arc simply weighted by the factor 1 — g(a) 
(integrated over a together with the other standard time- 
dependent terms). 

More specifically for r c = 300 Mpc and a£/r c -C 1, ac- 
cording to our analytic results, the contribution to the 
LG peculiar velocity by a galaxy at £ is weighted by 
~ 1.21. For a£/?* c 3> 1, the galaxy's contribution is 
weighted by ~ 0.56. (These weights are found by per- 
forming the integral over the scale factor a in eq. (|10[) 
and they are defined relative to the predicted GR con- 
tribution). In Figure [5] the analytic limit for a£/r c <C 1 
agrees well with the numerical results. We don't expect 
to reproduce the a£/r c ^> 1 limit since our sample starts 
at ~ 5 Mpc and the analytical results are not reliable for 
r c < 50 Mpc. 



FIG. 7: Components of the LG peculiar velocity as a function 
of size of sphere-of-influence under degravitation formulation. 
The vertical axis is normalized such that at the limiting dis- 
tance of the sample (140 Mpc), the magnitude of the GR-case 
velocity is unity. The graphs show the cartesian components 
of the LG peculiar velocity in galactic coordinates. The solid 
line represents the GR reconstruction; the dotted and dashed 
lines represent the peculiar velocities reconstructed for two 
different scale lengths: r c = 20 Mpc and 300 Mpc, respec- 
tively. While r c = 20 Mpc is not physically meaningful, it 
is insightful to understand how a theory of this form could 
change the LG peculiar velocity. Note that the vertical axes 
of the different cartesian components do not have the same 
scales. 



Figure O shows that the angular deviation from GR 
is bigger for smaller r c , with a value of ~ 6° at r c ~ 
20 Mpc. However, this small value of r c is inconsistent 
with the matter power spectrum and other cosmological 
observables [3^|; for values of r c in the relevant regime 
300—600 Mpc, the angular deviation is less than a degree; 
therefore these higher-dimensional models cannot explain 
the observed misalignment angle. The angular deviation 
goes to zero for a£/r c <C 1 since the direction of the LG 



12 



dipolc does not change when the contributions from all 
the masses are uniformly rescaled. Similarly, we expect 
the angular deviation to go to zero if the domain were 
extended to values of r c such that a£/r c ^> 1. 

Since the survey is limited to £ < 140 Mpc, it is 
straightforward to understand the numerical results for 
values of 300 Mpc < r c < 600 Mpc. As r c increases, mod- 
ifications to GR become relevant at later times, leading to 
the integrated time-dependent correction (1 — g(a)) hav- 
ing a smaller effect. In addition, larger r c causes more of 
the masses to be weighted by a common factor, leading 
to less rotation of the dipole. 

We can again use the components of the peculiar veloc- 
ity in concentric spheres, as shown in Figure [3 to under- 
stand how different masses contribute to the LG peculiar 
velocity. The detailed analysis of Figure|3]in the previous 
subsection, can be reapplied to Figure [7l When consid- 
ering values of r c close to the distance of the Virgo Clus- 
ter, the relative suppression of more distant groups shifts 
the LG peculiar velocity towards Virgo, and also slightly 
closer to the CMB-implied direction. For r c = 300 Mpc, 
all galaxy groups satisfy a£/r c <C 1 and so the LG pecu- 
liar velocity is simply amplified by ~ 1.21 with negligible 
change in the direction. Note once again that since the 
amplitude of the change in angle is significantly less than 
the error on the reconstruction, there is little insight to 
be gained from the direction of the change. 



C. Summary 

As expected, the scale-dependent modifications to the 
Poisson's equation (j4j rotate the direction of the recon- 
structed LG peculiar velocity. However, for the broad 
set of models with realistic parameter choices that we 
have considered, the change in direction with respect to 
the standard GR result is of the order of a few degrees, 
which is significantly smaller than the ~ 40° needed to 
align the LG velocity with the CMB dipolc. 

From the analysis it is clear that, if there were a large 
undetected mass either inside or outside the survey vol- 
ume, which could rotate the dipolc in the right direction, 
modified gravity models would still predict a direction 
for the LG peculiar velocity within a few degrees of the 
GR prediction. In other words, modified gravity would 
not be significantly more efficient than GR in solving the 
angular discrepancy. 

From Figures [3] and [51 it is clear that for reasonable 
parameter choices, the magnitude of the LG peculiar ve- 
locity differs at most by a factor of 2 from the GR result. 
Since there are several uncertainties, like the bias, which 
are used to determine the magnitude of the LG peculiar 
velocity, we cannot discount modified gravity on the ba- 
sis of the magnitude of the reconstructed velocity alone. 
In other words, we can always find a parameter b of or- 
der unity that will allow the peculiar velocity of the LG 
to equal the magnitude of the peculiar velocity inferred 
from the CMB dipole. 



V. CONCLUSIONS 

The peculiar velocity of the Local Group, recon- 
structed from inhomogeneities in the density field, does 
not agree with the velocity inferred from the dipole 
anisotropy in the CMB. Using the HDC catalog of galaxy 
groups pj], EU constructed from 2MRS [HI, yj| , we in- 
vestigated, with both numerical and analytical methods, 
whether allowing for modifications of gravity on large 
scales, or interactions in the dark sector, can alleviate 
this discrepancy. Using an approximate expression for 
the growth factor, we have derived a general analytical 
formalism that relates peculiar velocities in real space to 
matter overdensities. This formalism applies to general- 
ized theories of gravity that reproduce GR at early times 
and predict a modified Poisson equation on sub-horizon 
scales (eq. 0J- Moreover, it reproduces numerical results 
for a broad range of parameter values. 

We have then specialized to two broad classes of mod- 
els, representing the main approaches to the phenomenon 
of cosmic acceleration, and reconstructed the correspond- 
ing LG peculiar velocity. The models we have chosen 
have the common feature of introducing scale-dependent 
modifications to the growth of perturbations, and there- 
fore have the potential of rotating the direction of the 
peculiar velocity. Specifically, we have considered scalar- 
tensor theories and models of dark energy coupled to dark 
matter as well as higher-dimensional gravity [19l . |20T ] . For 
reasonable values of each model's parameters, we find 
that the direction of the LG peculiar velocity changes by 
at most ~ 3° with respect to the standard GR result. 
This angle is much less than the required ~ 40° to co- 
incide with the direction of the LG dipole inferred from 
the CMB. 

The inability to reconcile the direction of the LG pecu- 
liar velocity with the CMB dipole, under either ACDM or 
modified gravity, continues to raise numerous questions. 
The selection-bias of magnitude-limited surveys prevents 
us from detecting all luminous matter inside the limit- 
ing survey volume; however, increasing the magnitude 
limit in 2MRS docs not significantly increase the galaxy 
count inside 7000 km s~ x ; moreover, since we considered 
galaxy groups, we expect to have included the most mas- 
sive structures inside 140 Mpc in our analysis. Further- 
more, while previous efforts have relied upon assump- 
tions about the masses of galaxies, we have used masses 
estimated dynamically from galaxy group members, and 
therefore have indirectly accounted for dark matter inte- 
rior to the observed galaxies in clusters. 

Even though 2MASS samples galaxies at very low 
galactic latitudes, it is possible that we have missed a 
massive contribution to the LG velocity located very close 
to the galactic plane. While an as-yet undetected con- 
tribution to the LG velocity (e.g., a poorly sampled, dis- 
tant, massive supercluster, or — as suggested by [46| - 
a nearby galaxy masked by the galactic plane) could rec- 
oncile the issue, we find that GR and modified gravity 
are expected to give directions for the LG velocity that 
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are within a few degrees of each other. A modification 
of the laws of interaction would therefore not play a key 
role if this was the case. 

On the other hand, the significant discrepancy between 
the two dipoles leads us to question the assumptions that 
are built into the reconstruction. In this and previous 
work, it is assumed that the luminous matter traces the 
dark matter. While we have attempted to estimate the 
cluster masses dynamically (using luminous matter as 
tracers), we can only estimate interior mass. Our as- 
sumption of linear bias implies that the estimated cluster 
mass is proportional to the dark matter mass, and that 
this constant of proportionality does not depend on the 
properties of the cluster, or on rcdshift, which may not 
be a valid assumption. 

We have shown that changing the scale-dependence 
of the growth of perturbations does alter the recon- 
structed direction of the LG dipole, although not in an 
amount sufficient to explain the discrepancy with the 
CMB. While we have assumed a linear bias between visi- 
ble and dark matter, a scale-dependent bias might have a 
similar effect to that of modified gravity discussed above, 
but subject to different observational constraints. It is 
evident that a scale-dependence in the bias parameter 
will alter the direction of the reconstructed LG velocity; 
however, from our analysis, we cannot conclude whether 
or not realistic models of scale-dependence in the bias 
could account for the angular discrepancy. 

An alternative explanation is a limit in our understand- 
ing of cosmology. In a ACDM universe, we expect that 
the motions of galaxies (and clusters of galaxies) merge 
with the Hubble flow at ~ 150 Mpc scales. However 
measurements of the power asymmetry in the CMB [47j 
and of bulk flows [H, [4^| , suggest an inhomogeneous dis- 
tribution of matter on much larger scales. This large- 
scale power, while inconsistent with ACDM, could ex- 



plain the discrepancy between dipoles. Indeed, the dipole 
reconstructed from local inhomogeneities assuming GR, 
in conjunction with a bulk flow in the direction suggested 
by both [H| and [49[, can result in a net flow in agree- 
ment, within the estimated errors, with the direction ob- 
served in the CMB (with an appropriate choice for the 
bias parameter). Further discussion on this subject, how- 
ever, is beyond the scope of this work (see [5(3 for a de- 
tailed analysis). 

Peculiar velocities probe the matter density on a wide 
range of scales, and thus present a unique test of gravity 
on large sub-horizon scales. While the effect of modified 
gravity on the LG peculiar velocity is minimal, study- 
ing the impact of modified gravity on the local flow-field 
could potentially act as a more distinguishing test of 
gravity. Further study on this subject would be highly 
insightful, but is, once again, beyond the scope of this 
work. 

While in this work we cannot deduce the correct expla- 
nation for the discrepancy between the peculiar velocity 
of the LG induced by inhomogeneities in the density field 
and the velocity expected from the dipole in the CMB, 
we have demonstrated that modified theories of gravity 
alone cannot provide a solution to this problem. 
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